Efficient Computation of Casimir Interactions between Arbitrary 3D Objects 
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We introduce an efficient technique for computing Casimir energies and forces between objects of 
arbitrarily complex 3D geometries. In contrast to other recently developed methods, our technique 
easily handles non-spheroidal, non-axisymmetric objects and objects with sharp corners. Using our 
new technique, we obtain the first predictions of Casimir interactions in a number of experimentally 
relevant geometries, including crossed cylinders and tetrahedral nanoparticles. 
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Since the dawn of the modern era of precision Casimir- 
force measurements some 10 years ago |T], Casimir forces 
have been measured in an increasingly wide variety of ex- 
perimental geometries, including plate-plate [2], sphere- 
plate sphere-comb and cylinder-cylinder [5] con- 
figurations. A recent experiment [B] finds evidence of 
Casimir interactions in a commercially fabricated MEMS 
device, and there is every reason to believe that accurate 
modeling of Casimir forces in complex geometries will be 
a critical ingredient in the design of future commercial 
MEMS technologies [3 |HJ [3] ■ 

Theoretical developments have largely failed to keep 
pace with this rapidly advancing experimental reality. 
Until recently, theoretical considerations were restricted 
to simple models in highly idealized geometries [10] , such 
as infinite parallel plates or infinite parallel cylinders. 
Techniques for predicting Casimir forces in more gen- 
eral geometries are clearly needed if theory is to confront 
the growing wealth of experimental Casimir data and the 
design challenges of future MEMS devices. 

Important recent progress on general methods has 
been made by a number of authors |TTJ HH [131 IHl 
[T5] . Rodriguez et. al. jTT] applied standard tech- 
niques of computational electromagnetics to develop a 
method capable of predicting Casimir forces between ar- 
bitrary two-dimensional objects, i.e. infinitely extended 
objects of arbitrary 2D cross section. Emig et. al. 
(EGJK) [12] demonstrated an efficient algorithm for pre- 
dicting Casimir energies of interaction between com- 
pact 3D objects of spheroidal or nearly-spheroidal shape. 
Both of these techniques have since been applied to pre- 
dictions of Casimir interactions in new geometries |16l 
fTZl [18] . To date, however, a practical computational 
scheme for predicting Casimir forces between arbitrary 
non-spheroidal three-dimensional objects has remained 
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elusive. 

In this work, inspired by the ideas of Refs. [111 and [T^ 
we introduce a new technique for calculating Casimir 
forces between objects of arbitrary 3D shapes, including 
efficient handling of non-spheroidal, non-axisymmetric 
objects and objects with sharp corners. As an immediate 
application of our method, we present the first predic- 
tions of Casimir interactions in a number of experimen- 
tanlly relevant geometries, including crossed cylindrical 
capsules [S] and tetrahedral nanoparticles [HU (25] . 




FIG. 1: Surfaces of compact conducting objects are dis- 
cretized |19) into planar triangles, which are used to construct 
localized vector-valued basis functions describing the surface 
current distribution. As shown in the inset, each basis func- 
tion describes currents fiowing from a vertex of one triangle 
to the opposite vertex of an adjacent triangle, and vanishes 
on all other triangles j20| . 
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The method. The starting point of the EGJK method 
is a path-integral expression [23j for the Casimir energy 
of perfectly-conducting compact objects: 
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2^ 



dK log 
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Z{k)= J pj(x) (2) 

where the functional integration extends over all possible 
current distributions J(x) on the surfaces of the objects 

and where — [l+^V® V] 4^|x_x'| dyadic 
Green's function at (Wick-rotated) frequency —iui ~ en. 
Zoo in (jlj is -2 computed with all objects removed to 
infinite separation. 

Following EGJK, to obtain a tractable expression for 
Z for a given collection of No objects we now proceed 
to expand the current distribution in a discrete set of 
basis functions, J(x) = ^ Jimfim(x), where i = 1- ■ ■ Ng 
ranges over the objects in the geometry and to = 1 • • • A'^^; 
runs over a set of Ni expansion functions defined for the 
ith object. With this expansion we have 

= (2^)"/2^- [detM(K)]-^/' (4) 
and ([T| becomes 

he p detM(Ac) 
27r 7o detMoo(K) 

where we have evaluated the multiple Gaussian integral 
in (|3| using standard path-integration techniques [2^ . 
Here the elements of the matrix M(k) are the interac- 
tions between the basis functions, Mq,^(k) = //fa • Qk. • 
dxdx', and in ([s]) is the (constant) Jacobian of the 
transformation I? J — * Y\ dJim, which cancels in the ratio 

EGJK took the expansion functions fi„i to be 
proportional to the spherical multipole moments 
{QE.irmQM,im} ^^'^ ^^^^ objcct. Although this 
choice leads to rapidly convergent and even analytically 
tractable series for spherical objects, it is not of practical 
use for general geometries, as the spherical multipole ex- 
pansion must be carried to high orders to represent source 
distributions on highly non-spherical objects. Ref. [25| 
for example, demonstrated poor convergence rates for 
spherical-harmonic representations of elongated objects 
and objects with corners or cusps. 

A more general strategy is to discretize the surfaces of 
the objects into planar triangles, as depicted in Figure 
[l] and to introduce localized basis functions in the spirit 
of finite-element and boundary-element methods. Fol- 
lowing standard practice in computational electromag- 
netics [20) . we choose basis functions defined on pairs of 



adjacent triangles. As indicated in the inset of Figure 1, 
to the mth internal edge in the surface discretization of 
object i we associate a localized basis function fim, which 
describes currents flowing on the two triangles that share 
that edge and vanishes on all other triangles. The matrix 
elements M^^ reduce to finite surface integrals over plane 
triangles, the evaluation of which is greatly facilitated by 
the elegant methods of Ref. |26| and the overall accuracy 
of the computation may be systematically improved, at 
greater computational expense, by reducing the size of 
the triangles in the surface mesh. A further virtue of this 
choice of basis functions is that it eliminates the unwieldy 
distinction between interior and exterior fields in |12j . 

Casimir Energy. To compute the energy of interaction 
of our objects we write 
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where {A„},{A^} are the eigenvalues of M,Moo- 

Casimir Force. The z-directed force on the ith object 
in our geometry is obtained by differentiating ([5| with 
respect to displacement of that object: 

he f°° , d , detM(K) 
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The integrand may be conveniently evaluated using the 
identity g7] 
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where {a„} are the eigenvalues of the generalized eigen- 
value problem 



dm 

dz. 
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The n integrals in ([S]) and ([7| may be evaluated using 
standard numerical quadrature techniques. 

The dimension of the matrix M is the number of basis 
functions in our expansion of the surface current distri- 
bution, which in turn is proportional to the number of 
triangles in our surface discretization (see Figure For 
matrices of moderate dimension, N < 5000, correspond- 
ing to moderately fine mesh discretizations (sufficient for 
accurate treatment of all geometries considered here) , the 
eigenvalues in (|6| and ([8| may be computed using direct 
methods in 0{N^) time. For finer meshes it may be pos- 
sible to use iterative eigenvalue solvers and matrix spar- 
sification techniques [251 HSl [3D] to reduce the complexity 
to 0{N In N); this will be discussed in a subsequent pub- 
lication. 

Results. We now apply our technique to the prediction 
of Casimir interactions in a number of 3D geometries. 

Casimir energy of identical spheres, parallel capsules, 
and crossed capsules. Figure 2 shows plots of Casimir 



3 




1e-05 I ' 1 1 ' 1 1 ' ^ 

0.5 1 1.5 2 2.5 3 3.5 4 

Z/R (surface-surface separation) 



FIG. 2: Casimir energy vs. separation distance for three pairs 
of perfectly conducting objects: identical spheres (red curve), 
identical capsules with parallel axes (blue curve), and iden- 
tical capsules with perpendicular axes (green curve). The 
hollow circles represent the sphere-sphere data of Ref. 1121 

01 I 1 1 1 1 1 i 1 i 1 1 




0001 I 1 1 1 1 ' ' 1 ^ ' 1 

3 4 5 6 7 8 9 10 11 12 13 

L/R 

FIG. 3: Magnitude of attractive Casimir force between 
crossed capsules of radius 7? as a function of capsule length, for 
surface-surface separations Z = 2R and Z — 4R. The dashed 
lines denote the asymptotic (L oo) limits of the force, as 
roughly extrapolated from our finite-L data. The solid lines 
are fits of the large-L data to the form a + b/L, although 
our data are insufficient to establish the precise asymptotic 
L-dependence of the force. 

energy versus separation distance for three pairs of per- 
fectly conducting objects: identical spheres of radius R, 
identical capsules with parallel axes, and identical cap- 
sules with perpendicular axes ("crossed capsules.") (A 
"capsule" is a cylinder of radius R with hemispherical 
endcaps, of total length L = 6R in this case.) The en- 
ergy curve for the crossed capsules interpolates between 
the curve for the spheres at short distances and the curve 
for the parallel capsules at large distances. We may un- 
derstand this result as follows: At short distances, the 



interaction is dominated by contributions from the por- 
tions of the surfaces that lie in closest proximity to each 
other. For the crossed capsules, this region of nearest 
proximity is restricted to a length ~ R around the cen- 
ters of the capsules and hence exhibits the same scaling as 
the region of nearest proximity for the two spheres in the 
short-distance limit. In contrast, for parallel capsules, 
the strongly-interacting region extends down the entire 
length of the capsule and hence scales with an additional 
factor of L at short distances. At large distances, the 
interaction becomes pointlike and the energy of interac- 
tion between identical capsules becomes independent of 
the orientation. 

Casimir Force vs. Length For Crossed Capsules. The 
Casimir force between crossed cylinders, in the limit in 
which the length L of the cylinders is much longer than 
their surface-surface separation Z, is experimentally ac- 
cessible [5] and interesting in its own right as a finite 
interaction between effectively infinite objects. Figure 
3 plots the magnitude of the attractive Casimir force 
between crossed metallic capsules of radius R, at fixed 
surface-surface separations oi Z — 2R and Z = 4R, as 
a function of the capsule length L. In the limit L —^ oo 
our capsules become infinite cylinders and the force ap- 
proaches a ^-dependent constant (denoted by the dashed 
lines in the figure.) 

An interesting question is how rapidly the force ap- 
proaches the infinite-cylinder limit as L — > oo at fixed Z. 
A numerical determination of the precise asymptotic in- 
dependence of the force requires an iterative fast-solver 
version of our method that is capable of handling ma- 
trices M of larger dimension; this extension will be dis- 
cussed in a future publication. 

Tetrahedral nanoparticles. Several groups have suc- 
ceeded in producing tetrahedral nanoparticles of dimen- 
sions ~ 10-50 nm pT| 122], and Casimir forces may dom- 
inate the interactions between electrostatically neutral 
particles of this type. Our method allows the first pre- 
dictions of Casimir interactions between perfectly con- 
ducting tetrahedral particles. (In a real system with par- 
ticles of these sizes, the finite conductivity of the met- 
als cannot be neglected. We are currently implementing 
an extension of our method to treat objects of arbitrary 
permittivities. Even for the idealization of perfectly con- 
ducting nanoparticles, however, it is interesting to ana- 
lyze the orientation dependence of the Casimir energy, 
as well as to demonstrate the practical ability of our 
method to treat geometries with sharp corners and non- 
axisymmetric shapes.) Figure |4] depicts a contour plot of 
Casimir energy vs. orientation angles for two perfectly 
conducting tetrahedra originally separated by a distance 
D — 2L in the y direction, where L is the tetrahedron 
edge length. To one of the tetrahedra we apply a rotation 
through an angle ip about the z axis followed by a rota- 
tion through an angle 6 about the y axis; the fixed origin 
of these rotations, depicted as the origin of the coordi- 
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FIG. 4: Contour plot of Casimir energy vs. orientation angles for tetrahedral nanoparticles separated by a distance D = 2L. 



nate axes in Figure |4] is the point lying a distance H/2 
below the apex of the tetrahedron, where H ~ \/3L/2 is 
the height of the tetrahedron. The contour plot reveals 
a clear minimum at {9, (f))=(0,60°), corresponding to the 
closest approach of a vertex of the base of the rotated 
tetrahedron to the unrotated tetrahedron, as might be 
expected for an attractive interaction. 

Conclusion. In conclusion, we have developed a 
general method for computing Casimir energies and 
forces between objects of arbitrary three-dimensional 
shapes, enabling efficient treatment of non-spheroidal, 
non-axisymmetric objects and objects with sharp cor- 
ners. Using our method, we have predicted Casimir inter- 
actions in a variety of experimentally relevant geometries 
that would be challenging to handle by previous meth- 
ods. 

Acknowledgements. We are grateful to T. Emig for 
providing the raw data from Ref. [121 The authors ac- 
knowledge the support of the Interconnect Focus Cen- 
ter, one of five research centers funded under the Focus 
Center Research Program, a DARPA and Semiconduc- 
tor Research Corporation program. In addition, the au- 
thors are grateful for support from the Singpore-MIT Al- 
liance Computational Engineering flagship research pro- 
gram. This work was supported in part by a Depart- 
ment of Energy (DOE) Computational Science Fellow- 
ship under grant DE-FG02-97ER25308 (ARW), by the 
Army Research Office through the ISN under Contract 
No. W911NF-07-D-0004, and by the MIT Ferry Fund. 



[1] S. K. Lamoreaux, Phys. Rev. Lett. 78, 5 (1997). 

[2] G. Bressi, G. Carugno, R. Onofrio, and G. Ruoso, Phys. 

Rev. Lett. 88, 041804 (2002). 
[3] D. E. Krause, R. S. Decca, D. Lopez, and E. Fischbach, 

Phys. Rev. Lett. 98, 050403 (2007). 
[4] H. B. Chan et al., Phys. Rev. Lett. 101, 030401 (2008). 
[5] T. Ederth, Phys. Rev. A 62, 062104 (2000). 
[6] R. Ardito et al., IEEE Sensors 2008 Conference p. 90 



[7 

[8 
[9 

[10 
[11 

[12 

[13 
[14 

[15 

[16 
[17 

[18 

[19 

[20 

[21 

[22 

[23' 
[24 

[25 

[26 

[27 
[28 



(2008). 

F. Capasso, J. N. Munday, D. lannuzzi, and H. B. 
Chan, IEEE. J. Select. Topic Quantum Electron. 13, 400 
(2007). 

J.-G. Guo and Y.-P. Zhao, J. Microelectromechanical 
Systems 13, 1027 (2004). 

E. Oilier et al., IEEE. International Conference 
on Nano/Micro Engineering and Molecular Systems, 
Bankkok, Thailand, 16-19 Jan. 2007 p. 180 (2007). 
H. B. G. Casimir, Koninkl. Ned. Adak. Wetenschap. 
Proc. 51, 793 (1948). 

A. W. Rodriguez, M. Ibanescu, D. lannuzzi, J. D. 
Joannopoulos, and S. G. Johnson, Phys. Rev. A 76, 
032106 (2007). 

T. Emig, N. Graham, R. L. Jaffe, and M. Kardar, Phys. 
Rev. Lett. 99, 170403 (2007). 

S. Pasquah and A. C. Maggs (2008), arXiv:0801.4385. 
A. Lambrecht and V. N. Marachevsky, Phys. Rev. Lett. 
101, 160403 (2008). 

H. Gies and K. Khngmiiller, Phys. Rev. Lett. 97, 220405 
(2006). 

A. Rodriguez et al., Phys. Rev. Lett. 101, 190404 (2008). 
A. W. Rodriguez, J. D. Joannopoulos, and S. G. Johnson, 
Phys. Rev. A 77, 062107 (2008). 

T. Emig, N. Graham, R. L. Jaffe, and M. Kardar (2008), 
0811.1597. 

C. Geuzaine and J.-F. Remade, Int. J. Numer. Methods 
Eng. (2009), in press. 

S. M. Rao, D. R. Wilton, and A. W. Glisson, IEEE Trans. 

Antennas Propagat. AP-30, 409 (1982). 

T. S. Ahmadi, Z. L. Wang, T. C. Green, A. Henglein, 

and M. A. El-Sayed, Science 272, 1924 (1996). 

R. Jin, Y. Cao, C. A. Mirkin, K. L. Kelly, G. C. Schatz, 

and J. G. Zheng, Science 294, 1901 (2001). 

H. Li and M. Kardar, Phys. Rev. Lett. 67, 3275 (1991). 

S. Weinberg, The Quantum Theory of Fields, Vol. 1 

(Cambridge University Press, 1995), chapter 9. 

S.-H. Kuo, B. Tidor, and J. White, J. Emerg. Technol. 

Comput. Syst. 4, 1 (2008), ISSN 1550-4832. 

D. J. Taylor, IEEE Trans. Antennas Propagat. 51, 1630 
(2003). 

T. Emig, Europhys. Lett. 62, 466 (2003). 
Z. Zhu, B. Song, and J. White, in Proceedings of the De- 
sign Automation Conference (Anaheim, CA, 2003), pp. 
712 -717. 



5 



[29] V. Rokhlin, Journal of Computational Physics 60, 187 

(1985). 

[30] W. Hackbusch and Z. P. Nowak, Numerische Mathematik 
54, 463 (1989). 



